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Abstract 

Thiamethoxam has been used as a major insecticide to control the B-biotype sweetpotato whitefly, 
Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae). Due to its excessive use, a high level of 
resistance to thiamethoxam has developed worldwide over the past several years. To better 
understand the molecular mechanisms underlying this resistance in B. tabaci, gene profiles 
between the thiamethoxam-resistant and thiamethoxam-susceptible strains were investigated 
using the suppression subtractive hybridization (SSH) library approach. A total of 72 and 52 up- 
and down-regulated genes were obtained from the forward and reverse SSH libraries, 
respectively. These expressed sequence tags (ESTs) belong to several functional categories based 
on their gene ontology annotation. Some categories such as cell communication, response to 
abiotic stimulus, lipid particle, and nuclear envelope were identified only in the forward library of 
thiamethoxam-resistant strains. In contrast, categories such as behavior, cell proliferation, 
nutrient reservoir activity, sequence-specific DNA binding transcription factor activity, and 
signal transducer activity were identified solely in the reverse library. 

To study the validity of the SSH method, 16 differentially expressed genes from both forward 
and reverse SSH libraries were selected randomly for further analyses using quantitative real- 
time PCR (qRT-PCR). The qRT-PCR results were fairly consistent with the SSH results; 
however, only 50% of the genes showed significantly different expression profiles between the 
thiamethoxam-resistant and thiamethoxam-susceptible whiteflies. Among these genes, a putative 
NAD-dependent methanol dehydrogenase was substantially over-expressed in the thiamethoxam- 
resistant adults compared to their susceptible counterparts. The distributed profiles show that it 
was highly expressed during the egg stage, and was most abundant in the abdomen of adult 
females. 
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The sweet potato whitefly, Bemisia tabaci 
(Gennadius) (Hemiptera: Aleyrodidae), is one 
of the most widely distributed agricultural 
pests worldwide (P erring 2001), attacking 
agronomic, horticultural, and ornamental 
crops in subtropical and tropical agriculture, 
as well as in greenhouse production systems 
(Oliveira et al. 2001). It was first identified in 
China in the mid-1990s and then spread into 
more than 20 provinces within a very short 
time frame (Luo et al. 2002; Chu et al. 2005; 
Zhang et al. 2005; Chu et al. 2006). The 
phloem-feeding whitefly has caused severe 
crop losses through direct feeding, excretion 
of honeydew (which favors sooty mold 
development), and transmission of plant 
viruses (Jones 2003). 

Due to its severe damages, B. tabaci has been 
controlled predominantly with chemical 
insecticides. However, as a result of extensive 
application of synthetic insecticides, B. tabaci 
has developed a high degree of resistance to a 
wide range of insecticides, including 
carbamates, organophosphates, pyrethroids, 
insect growth regulators (IGRs), and 
neonicotinoids (Horowitz et al. 1988; 
Prabhaker et al. 1988; Horowitz et al. 1999; 
Nauen et al. 2002; Ahmad et al. 2002; Kranthi 
et al. 2002; Ma et al. 2007; Erdogan 2008; 
Roditakis 2009; Wang et al. 2010). 
Neonicotinoid insecticides are generally 
considered systemic and have excellent 
efficacy, long-lasting residual activity, and 
favorable safety profiles. For example, 
thiamethoxam, discovered and developed by 
the Novartis Crop Protection 
( www.novartis.com ), has played a crucial role 
in controlling B. tabaci and many other sap- 
sucking insect pests in China since its 
introduction in 2000. A high level of 
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resistance to thiamethoxam (100- and 900- 
fold), however, has already been reported in 
B- and Q-biotype of B. tabaci strains from 
Israel and Spain, respectively (Rauch and 
Nauen 2003; Horowitz et al. 2004). In China, 
both biotypes have developed high levels of 
resistance to imidacloprid and thiamethoxam 
in the field (Wang et al. 2010). 

In general, the safety and effectiveness of 
neonicotinoids have been attributed to their 
high affinity to nicotinic acetylcholine 
receptors (nAChRs). Consequently, resistance 
to neonicotinoids initially focused on the 
mutations in nAChRs (Liu et al. 2009). In 
addition, recent studies revealed that 
resistance of neonicotinoids in B. tabaci could 
be associated with an enhanced oxidative 
detoxification by cytochrome P450 
monooxygenases (Karunker et al. 2008; Wang 
et al. 2009). Through biochemical 
characterization of B-biotype thiamethoxam- 
resistant strains, cytochrome P450 
monooxygenase and carboxylesterase were 
found to be responsible for the thiamethoxam 
resistance in whiteflies (Feng et al. 2008, 
2010). 

Suppression subtractive hybridization (SSH) 
is a RNA-based method for identifying genes 
with unknown function, especially in species 
that lack primary genomic resources 
(Diatchenko et al. 1996; Lu and Wan 2008). 
This method has already been used to better 
understand the genetic basis of insecticide 
resistance, such as Aedes aegypti resistance to 
deltamethrin (Lertkiatmongkol et al. 2010) 
and Nilaparvata lugens resistance to 
triazophos (Bao et al. 2010). The SSH method 
has been applied to identify genes related to 
viral infection (Li et al. 2011) and heat-shock 
(Lii and Wang 2008) in B. tabaci as well. In 
this study, gene expression profiles between 
the thiamethoxam-resistant and 
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thiamethoxam-susceptible B. tabaci were 
investigated by both SSH and qRT-PCR 
analyses. Combined results give us a unique 
perspective in regards to the development of 
neonicotine resistance in the B-biotype B. 
tabaci. 

Materials and Methods 

Bemisia tabaci strains 

The B-biotype B. tabaci susceptible strain 
(TH-S) and resistant strain (TH-R) were the 
same populations as described previously 
(Feng et al. 2008, 2010). Before sample 
collection, a leaf-dip bioassay (Feng et al. 
2008) was conducted to confirm that the 
resistance factor (LC 50 (TH-R)/LC 50 (TH-S)) 
was at least 70-fold. About 3000 adult 
whiteflies from TH-R were treated with 2000 
mg/L thiamethoxam (-LCso) to eliminate the 
heterozygous individuals. Then, the survivors 
were collected after 48 hours and designated 
as the TH-2000. A total of 300 TH-S and TH- 
2000 adults, respectively, were collected, snap 
frozen in liquid nitrogen for three hours, and 
transferred to a -80 °C freezer for long-term 
storage. Different developmental stages, such 
as eggs, 3 rd instar larvae, and two-day-old 
unmated adult females, and various tissues 
including head, thorax, abdomen, and wing of 
a two-day-old unmated adult female were 
collected to study the distribution profiles of 
genes of interest. 

Total RNA isolation and reverse 
transcription 

Total RNAs from both TH-2000 and TH-S 
adults were extracted using Trizol (Invitrogen, 
www.invitrogen.com ) following manufacturer 
protocol. The resulting total RNA was re- 
suspended in nuclease-free water, and 
quantified by the Nanodrop 2000 (Thermo 
Scientific, www . thermoscientific .com ). The 
first-strand cDNA and ds-cDNA were 



synthesized using SMARTer™ PCR cDNA 
Synthesis Kit (Clontech, www.clontech.com ) 
and later, the ds-cDNA were purified with 
QIAquick PCR Purification Kit (QIAGEN, 
www.quiagen.com ). 

Construction of the SSH cDNA library 

The SSH procedure was carried out using a 
PCR-Select™ cDNA Subtraction Kit 
(Clontech) following manufacturer protocol. 
A forward SSH library was constructed to 
isolate the up-regulated genes of the TH-2000 
whitefly strain. The forward SSH library was 
used to identify clones in which the TH-2000 
cDNA was used as the tester and the TH-S 
cDNA as the driver. In addition, a reverse 
SSH library was constructed to detect the 
down-regulated genes of the TH-2000 
whitefly strain. The reverse SSH library was 
used to identify clones in which the TH-S 
cDNA was used as the tester and the cDNA 
from TH-2000 as the driver. After 
hybridization, the subtracted cDNA were 
ligated into the pGEM-T vector (Promega, 
www.promega.com ) and transformed into 
Escherichia coli competent cells through 
electroporation. 

DNA sequencing and EST analysis 

Positive clones were selected by conventional 
blue-white screening. White clones were 
randomly selected from both forward and 
reverse libraries. The positive clones were 
further validated by colony PCR using nested 
PCR primers provided in the kit. The resulting 
products were subjected to the direct 
sequencing with Ml 3 primers. The vector 
sequences were removed through a Perl script 
and checked through VecScreen 
( http://www.ncbi.nlm.nih.gov/VecScreen/Vec 
Screen.html ). Then, the remaining high- 
quality EST sequences were analyzed in the 
GenBank non-redundant (nr) database with 
BLASTX. A sequence was considered as 
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significantly matched when the E-value was < 
10~ 5 . Functional annotation was carried out in 
the Swiss-Prot 
( http://expasy.org/people/swissprot.html) . 
Gene Ontology (GO) terms were extracted 
and counted using map2slim and Perl scripts. 

Quantitative real time PCR 

Up to 150 (three biological replicates, n = 50) 
TH-S and TH-2000 adults, respectively, were 
collected for the qRT-PCR analysis. 
Approximately 0.5 ug of total RNA was used 
as a template to synthesize the first-strand 
cDNA using PrimerScript RT reagent Kit 
(Takara Bio Inc., www.takarabio.com ) 
following manufacturer protocol. The 
resulting cDNA was diluted to a working 
concentration of O.lug/uL for the subsequent 
qRT-PCR analysis. To validate the 
differentially expressed genes detected by the 
SSH approach, 16 expressed sequence tags 
(ESTs) representing 11 putatively up- 
regulated genes and five down-regulated 
genes were randomly selected. The qRT-PCR 
primers were designed using Primer3 



( http://frodo.wi.mit.edu/primer3 ) (Table 1). 
The cycling parameters were as follows: 95 
°C for 30 sec, followed by 40 cycles of 95 °C 
for five sec and 62 °C for 34 sec, and ended 
with a melting curve analysis (65 °C to 95 °C 
in increments of 0.5 °C every five sec) to 
check for nonspecific product amplification. 
Relative gene expression was calculated using 
the 2" AACt method, p-actin was used as the 
internal reference gene to eliminate sample-to- 
sample variations in the initial cDNA samples. 

Results 

Differential screening and EST sequencing 

Based on the results of the differential 
screening, all 507 cDNA clones were 
randomly picked and sequenced from these 
two libraries. Specifically, 298 clones were 
from the forward library representing up- 
regulated genes, and 209 clones were from the 
reverse library representing down-regulated 
genes. After trimming, 127 and 63 high- 
quality ESTs from the forward and reverse 
library, respectively, were obtained. 
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Figure I. Gene Ontology (GO) classification of the differentially expressed expression ESTs in the forward and reverse SSH 
libraries. Based on the GO terms, the ESTs are categorized into putative functional groups. The black and white columns 
represent the up- and down-regulated genes from forward and reverse libraries, respectively. High quality figures are available 
online. 
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Figure 2. Quantitative real-time PCR validation. The gene expression level of 1 6 randomly selected ESTs, including I I from 
SSH forward library (the black column) and five from SSH reverse library (the white column), was tested using qRT-PCR. The 
relative gene expression level in the resistant strains was normalized to the susceptible TH-S whiteflies. Data are presented as 
mean ± SE. Asterisks denote significant gene expression differences between the resistant and susceptible whiteflies, as 
determined by the pairwise t-tests (* p < 0.05, ** p < 0.0 1, LSD t-test). High quality figures are available online. 



Functional annotation 

The BLASTX results showed that among the 
127 clones from the SSH forward library, 72 
ESTs (56.7%) had significant matches (E- 
value < 10" 5 ) to known or predicted genes in 
GenBank, and these clones could be 
assembled into 61 distinct sequences 
(Supplementary Table 1). For the 63 clones 
from the SSH reverse library, 52 ESTs (83.8 
%) had significant matches (E -value < 10" 5 ) 
with the database, among which 39 distinct 
sequences were identified (Table S2). 

Based on the Gene Ontology (GO) terms, 
these distinct sequences were functionally 
annotated (Figure 1). Majority of the 
biological processes (such as cell cycle, 
metabolic process, response to external 
stimulus, response to stress, and transport), 
some cellular components (such as cell, 
plasma membrane, and protein complex), and 
some molecular functions (such as binding, 



catalytic activity, structural molecule activity, 
transcription regulator activity, and transporter 
activity) were presented in both libraries. 
Other GO terms, such as cell communication, 
response to abiotic stimulus, lipid particle, and 
nuclear envelope were identified only in the 
forward library. Vice versa, some GO terms 
such as behavior, cell proliferation, nutrient 
reservoir activity, sequence-specific DNA 
binding transcription factor activity, and 
signal transducer activity were identified only 
in the reverse library. 

qRT-PCR validation 

The qRT-PCR results from the randomly 
selected differentially expressed transcripts 
were, for the most part, consistent with the 
SSH results (Figure 2, Supplementary Tables 
1 and 2). For the 1 1 up-regulated genes, 9 of 
them were over-expressed in the resistant B. 
tabaci, whereas 4 out of 5 down-regulated 
genes were under-expressed. However, only 
50% of the genes exhibited significantly 
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Head Thorax Abdomen Wing 

Figure 3. Gene expression profiles of a putative NAD- 
dependent methanol dehydrogenase in Bemisia tabaci. (A) 
Relative gene expression between resistant (TH-2000) and 
susceptible (TH-S) Bemisia tabaci adults. (B) Relative gene 
expression at different developmental stages including egg, 
3 rd instar larvae, and two-day-old unmated adult females. 
(C) Relative gene expression among different tissues 
including head, thorax, abdomen, and wing of two-day-old 
unmated adult females. High quality figures are available 
online. 



different expression profiles between resistant 
and susceptible whiteflies (Figure 2). Most 
notably, F-THSS58, a putative NAD- 
dependent methanol dehydrogenase EST, was 
over-expressed ~ 12-fold in the resistant TH- 
2000 whiteflies in comparison to the 
susceptible TH-S strains (Figure 2, 
Supplementary Table 1). To characterize this 
newly identified B. tabaci dehydrogenase 
gene, its expression profiles at different 
developmental stages and different tissues 
were examined (Figure 3). In general, the 
transcript level of this gene was much higher 
in the resistant whiteflies (Figure 3A), was 
most abundant in egg stage (Figure 3B), and 
was much higher in the abdomen of adult 
female than in any other tissues (Figure 3C). 
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Discussion 

Thiamethoxam, a second-generation 
neonicotinoid insecticide (Maienfisch et al. 
2001), has been used extensively for the 
sustainable management of B. tabaci in 
horticultural and other cropping systems 
(Nauen and Denholm 2005). However, like 
many other neonicotinoid insecticides, B. 
tabaci has developed a high degree of 
resistance to thiamethoxam under the 
laboratory selection (Feng et al. 2008, 2010), 
as well as in the field (Elbert and Nauen 2000; 
Horowitz et al. 2004; Wang et al. 2010) in the 
past decade. The molecular mechanism 
governing the thiamethoxam resistance in B. 
tabaci, however, has yet to be fully 
understood. In this study, the molecular basis 
of thiamethoxam resistance in B. tabaci was 
investigated using the SSH cDNA library 
approach. About 72 and 52 differentially 
expressed ESTs were obtained from forward 
and reverse libraries, respectively, 
representing up-regulated and down-regulated 
genes. The differentially expressed genes 
between the thiamethoxam-resistant and - 
susceptible B. tabaci include, but not limit to, 
cell communication, response to abiotic 
stimulus, response to stress, lipid particle, 
nuclear envelope, cell proliferation, and 
nutrient reservoir activity, etc. The accuracy 
of the SSH method was partially confirmed by 
the qRT-PCR analysis, with only 50% of the 
randomly selected ESTs showed significant 
differences. Similar to microarray analysis, 
RNA-based SSH method tends to generate 
false positives. Consequently, results from 
both analyses need to be validated by qRT- 
PCR. 

Previous mechanistic studies suggested that 
neonicotinoid resistance could be associated 
with enhanced metabolic detoxification by 
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cytochrome P450 monooxygenases (Zhao et 
al. 2000; Rauch and Nauen 2003; Honda et al. 
2006; Karunker et al. 2008; Wang et al. 
2009). In this study, only one P450-like EST 
(F-THSS19), that has the highest similarity 
with a tobacco cytochrome P-450-like gene 
(Accession No. BAA10929.1), was 
significantly over-expressed in the resistant 
whiteflies (1.50 ± 0.05, p<0.01, Table 2). Due 
to the limited resolution and coverage of this 
SSH method, it is not uncommon that some of 
the genes potentially involved in the 
thiamethoxam resistance in whiteflies were 
not included. It is worth noting, however, that 
a NAD-dependent methanol dehydrogenase- 
like EST from B.tabaci was substantially 
overexpressed in the resistant whiteflies 
(11.56 ±0.57, pO.Ol, Table 2). 

Dehydrogenases including farnesol 
dehydrogenase, succinic semialdehyde 
dehydrogenase (SSADH), aldehyde 
dehydrogenase, glutamate dehydrogenase, and 
methanol dehydrogenase can oxidize a 
substrate by transferring one or more hydrides 
(H~) to an acceptor, usually NAD + /NADP + or 
a flavin coenzyme such as FAD or FMN. 
They are involved in various physiological 
and biochemical processes. In mammals, 
SSADH is thought to be responsible for the 
degradation of the inhibitory neurotransmitter 
GABA in the central nervous system (Blaner 
and Churchich 1979; Chambliss et al. 1995). 
SSADH homologues have been cloned and 
expressed in the parasitic insects Lucilia 
cuprina and Ctenocephalides felis (Rothacker 
et al. 2008). In addition, NADP + -dependent 
farnesol dehydrogenase was found to be 
involved in the juvenile hormone synthesis in 
mosquito (Mayoral et al. 2009). The NAD- 
dependent methanol dehydrogenase found in 
this study shed new light on the molecular 
understanding of thiamethoxam resistance in 
whiteflies. Based on these results, future 
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studies involving cloning and functional 
characterization of this NAD-dependent 
methanol dehydrogenase are warranted to 
elucidate its role in the whitefly thiamethoxam 
resistance. 
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Table I. Primers used for the quantitative real-time PCR analysis. 



Putative Gene 


Primers 5'-3'* 


Temp (°C)** 


Product (bp) 


Alcohol dehydrogenase 


F-CAATGCGAGCTGCTCTGGA 


62 


102 


R (it 1 GUAAI GGACGAG IGGAA 




F-CATCAGTAATCGCGGATCGC 


62 


101 


R-CCATCCGCACTTTTTGCAAC 


(.ytoenrome 1 h jv— iiKe i isi 


F-TCTGCCCAGTGCTCTGAATGT 


62 


101 


R-TGACGAGGCATTTGGCTACC 


Heat shock protein 90 


F-CCTCCGGTTTTGCTCTTGAAG 


62 


104 


RTTCAACCATGACTGGCTCGTC 


Zinc finger protein 585A 


F-AGCCATCTGCCTGCTGATCTC 


62 


104 


R-CC 1AGATGGA1 1CGCCTGIGG 


Kinesin heavy chain 


F-TGATCGGCTCTTTGAGCAATG 


oz 


1 si 

151 


R-AATTCGTAGTCGCGGCGAG 


Alanine aminotransferase 


F-ACACAGCCTCATTGGCTGAAT 


oz 




R-CAGGAATATCCAAGCCCCAGT 


Amino acid transporter 


F- GGCATCACACTTCCGATTGC 


oz 




R-CCCAACATTGCTGAACGCA 


NAD-dependent methanol dehydrogenase 


F-TCGAGTCCri'CACAGCAlTGC 


62 


107 


R-GGCCCAAATAICGGAGAI GG 


Pyridoxal kinase 


i- - I GAGACGGAAG1 AGTGCCC'AC 


62 


116 


R-TCCTTCCGTTTCATGCTGTTC 


Ubiquitin-conjugating enzyme E2 


F-AGTCGGCGAACTTAACACGCT 


62 


102 


R-CGCCAGCCTTTCCTTTCTTC 


Eukaryotic translation initiation factor 3 subunit 7 


F-AAGCCCCTGTCAAAAATCGTG 


62 


136 


R-CCACAGCAGGTAAGGTCTGCA 


Vltellogenin-3 


F-AAAAATACGCCCGC'AACGA 


62 


144 


R-GGATGGAAGGTCCGTTTTGG 


CREB isofonn 2 


!- CTGAI 1GCIGGI I'GGCI I'GAA 


62 


117 


R-GGATGGTGGAGGAAAACAGGA 


Reverse transcriptase 


F-CGAGCCGACCAACTAAACCA 


62 


101 


R-ACACCTAGTGTTGCAGCCGG 


Vasodilator-stimulated phosphoprotein 


F-GGATGGTGATGGTGTGTCAGG 


62 


129 


R-CCGCACCCAACAAGCTCTT 


P-actin 


F- ACCGCAAGATTCCATACCC 


60 


129 


R-CGCTGCCTCCACCTCATT 



* F, forward primer; R, reverse primer; ** T, Annealing temperature. 
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Supplementary Table I. Up-regulated genes identified in the SSH forward library. 



Clone ID* 


Length (bpl 


Sequence Description 


Accession >'o. 


Hit Species 


E-valuc** 


Copv No. 


F-T1I SS 1 


663 


D-tvrosyl-tRNA(Tvr) deacvlase 1 


EFN74087.1 


Camponotus (loridunus 


1.00E-46 




F-TH SS 2 


311 


putative ribosomal protein L8e 


ABM55545.1 


Macunellicoccus hirsutus 


8. OOF. 33 


■ 


F-TH SS 3 


350 


non-ATPase reeulatorv subunit 


EFN85092.1 


Harpesnathos solicitor 


1.00E-07 




F-TH SS 4 


445 


6US nbosomal protein L7 


EFN75670.1 


Haroe&ncithos solicitor 


8 O0F-*>2 


, 


F-TH SS 5 


477 


60s nbosomal protein L9 


ADD20 176.1 


(ilassina runt si tans 


i 001-33 


j 


F-TH SS 6 


640 


AGAP002569-PA 


XP 312367.4 


Anopheles sambiae sir. PEST 


6.00E-25 


I 


F-TH SS 7 


817 


AGAP004890-PB 


XP 314288.4 


Anopheles gambiac str, PEST 


1.00E-138 




F-TH SS 8 


559 


Alcohol dehydrogenase 


XP 002429805.1 


Pediculus humanus corporis 


6.00E-71 


I 


F-TH SS 10 


379 


A1P synthase I 11 subunit 6 


YP 0868(15. 1 


Ih'rnisia tuhui i 


3.00K-34 


7 


F-TH SS 11 


670 


cathepsin B-912 


NP 001119612.2 


Acyrthosiphon pisum 


5.00E-30 


J 


F-TH SS 15 


561 




1K.R1 




8.00E-40 




F-TH SS 18 


306 


class III chitinase 


ACM45715.1 


Pvms pvrifoliu 


4.00E-27 




F-TH SS 19 


i26 


cytochrome P450 like THP 


nAA10929.I 


Nicotiana tahacum 


2.001.-2:? 


| 


F-TH SS 21 


703 


D-tyrosvl-tRNA(Tyr) deacylase 1 


EFN74087.1 


Camponotus /loridunus 


4.00E-47 




F-TH SS 22 


578 


tbx!7 


XP 002428011.1 


Pediculus humanus corporis 


6.00E-5 1 


| 


F-TH SS 23 


478 


tbnnin 


XP 002432855.1 


Pediculus humanus corporis 


5.00E-48 




I - I I 1 SS 26 


972 


GI13422 


XP 002008326.1 


Dmsophila mojavensis 


5.00F.-15 


| 


F-TH SS 27 


417 


GK20166 


XP 002060761.) 1 


Drvsophila willistoni 


2. 0011-30 




F-TH SS 28 


762 


Glycerophosphodiester phosphodiesterase 1 


EFN63344.1 


Camponotus floridanus 


6. OOF. -47 


| 


F-TH SS 29 


871 


OTP binding ADP ribosvlation factor like protein 1 protein 


EFN82659.1 


Harvevnalhos saltator 


7.00E-14 




K- TH SS 30 


461 


heat shock protein 90 


ACH85198 1 


tiemisia tahaci 


2 00F-V/ 




F-TH SS 31 


466 


TcasGA2 TC008350 


EFA02630.1 


Tribolium casturieum 


V00E 10 




F-TH SS 32 


384 


TcasGA2 TC009545 


EFA06626.1 


Tribolium castaneum 


5.00E-47 




F-TH SS 33 


346 


Importin subunit beta- 1 


EFN8 1834.1 


Harpegnathos saltator 


1.00E 30 




F-TH SS 34 


624 


Integrator complex subunit 3 


EFN89783.1 


Harpcv.nathos saltator 


1.00E-85 




F-TH SS 35 


396 


fCinesin heavy chain 


EFN7 1894.1 


Curnpunutus fluriilunus 


6. 00E -3 7 




F-TH SS 36 


757 


ribosonial protein, SI 1 


XP 001650870.1 


Aedes ueavpti 


7.00E-45 


J 


F-TH SS 37 


716 


Mitochondrial-processine peptidase subunit beta 


EFN65875.1 


Camponotus floridanus 


7. OOF 87 


J 


F-TH SS 38 


514 


multisvnthetase complex auxiliary component p38 


XP 00242746"' 1 


Pediculus humanus corporis 


5 OOF OX 


| 


T-Tll SS 19 


661 


Negative elongation factor A 


XP 002429474.1 


Pediculus humanus corporis 


5.00E-24 


] 


F-TH SS 40 


695 


Phosphoacelvlglucosaiiune tnutase 


EFN67555.1 


Camponotus /loridunus 


7.00E-17 




F-TH SS 41 


385 


acyl-Coenzyme A binding domain containing 4 


XP 002719487.1 


Oryctolagus cuniculus 


7.00E-06 


| 


F-TH SS 42 


767 


ADP-ribosvlation factor GTPasc activating protein 2-likc 


.XP 002740329.1 


Saccoslossus kowalevskii 


5 00E-09 


| 


F-TH SS 43 


47 1 


AGAP005404-PA 


XP 971530.1 


1 rihu/iu/n i ustiiiwum 


7 001 -1 i 


2 


F-TH SS 45 


923 


AGAP009657-PA 


XP 001947057.1 


Acyrthosiphon pisum 


5.00E-20 


j 


F-TH SS 46 


580 


AGAP0 10331 -PA 


XP 975279.2 


Tribolium castaneum 


5.00E-22 




F-TH SS 47 


599 


alanine aminotransferase 


XP 001948711.1 


Acvrthosiphon oisum 


1.00E-81 




F-TH SS 48 


458 


umino acid transporter 


XI' 01)1444972 1 


Acyrthosiphon pisum 


1 001-17 


j 


F-TH SS 50 


629 


aniinoacylase. putative 


XP 969212.1 


Tribolium castaneum 


3.00E-43 




F-TH SS 51 


716 


conserved hypothetical protein 


XP 001601024.1 


Nasonia vitripennis 


2.00E-L1 




F-TH SS 52 


298 


DEAD (Asp-GIu-Ala-Asp) box polypeptide 5 


XP 394723.2 


Apis mcllifcra 


2 00E-33 




l-TII ss 5:? 


750 


1)1' AI) hox A I P-dcpendcnt IiN A hclicase 


XP 001602K97.I 


Nasonia vitripennis 


7.00F-98 


■ 


F-TH SS 54 


376 


DEAD box ATP-dependent RNA helicase 


XP 974045.1 


Tribolium custuneum 


5.00E 43 




F-TH SS 55 


618 


GA15457-PA 


XP 001602440. 1 


Nasonia vitripennis 


6.00F-71 


j 


F-TH SS 56 


520 


GH10652O 


XP 001602045.1 


Nasonia vitripennis 


1.00E-52 




F-TH SS 57 


489 


MGC82386 protein 


XP 00lf>niK4K.I 


Nasonia vitripennis 


4.00F.-I4 




F-TH SS 58 


485 


NAD-dependent methanol dehydrogenase 


XP 001952147.1 


Acyrthosiphon pisum 


8. OOF- 60 




F-TH SS 59 


425 


nucleolar KKE/D repeat protein 


XP 001603746.1 


Nasonia vitripennis 


8.00F-44 




F-TH SS 60 


822 


pugilist CG4067-PA 


XP 623143.1 


Apis mellifera 


2.00E-80 




F-TH SS 61 


393 


Dvridoxal kinase 


XP 001607107.1 


Nasonia vitripennis 


6.00E-10 




F-TH SS 6? 


508 


ring finger protein 141 


XP 974067 1 


Tribolium castaneum 


9.00E-35 




F-TH SS 63 


599 


Type III alcohol dehydrogenase CG3425-PA 


XP 968236.1 


Tribolium castaneum 


6.00E-74 




F-TH SS 64 


544 


ubiquitin-coniugating enzyme E2 


XP 967918.1 


Tribolium castaneum 


1.00E-26 




F-TIi SS 65 


856 


Probable hvdroxvacid-oxoacid h-anshvdroecnasc. mitochondrial 


EFN79501.1 


Haroeenathos saltator 


4.00E-08 




F-TH SS 66 


395 


putative enolase 


ACT87785. 1 


Mezanola phvlla 


3 OOF- 10 




F-TH SS 67 


230 


REST corepressor 


EFN80005.1 


Harpeunathos saltator 


1.00E-14 




F-TH SS 68 


173 


ribosomal protein L8e 


ABF60235.1 


Leplinotarsa decemlineata 


4.00E-18 




F-TH SS 69 


456 


S-nhase kinase-associated protein 1 


ACOI1412.1 


Calieus roeercressevi 


4.00E-18 




F-TH SS 70 


548 


transposase IS630 


CA 1.4705 1.1 


Listonetta anguillarum semvar 02 


9 00E-06 




F-TH SS 71 


104 


transposase IS630 


CAL47051.1 


Listonella anauillarum serovar 02 


1.00E-06 




F-TH SS 72 


449 


Zinc finger protein 585A 


EFN63352.1 


Camponotus floridanus 


2.00E-54 


2 



* Genes randomly selected for the qRT-PCR validation study are in the shade; ** transcripts were putatively identified by 
homology from BLASTX search using a cutoff E-value less than IxlO- 5 . 
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Supplementary Table 2. Down-regulated genes identified in the SSH reverse library. 



Clone ID* 


Lenoth (bo) 


Sequence Description 


Accession No. 


HltSnedei 


E-value** 


Copy No. 


R-SS Til 1 


570 


MS ribosomal protein L27a 


XT 001892460.1 


Brueia malavi 


7.00E-54 


1 


R-SS TH 2 


466 


60S ribosomal protein L7 


EFN75670.I 


Harpesnalhos sallator 


5.00E-52 




R-SS 111 12 


375 


AT P synthase TO subunit 6 


YP 086805.1 


Bemisia (abaci 


3.00 L-J4 


T7I 

10 


n-ss m i.i 


902 


ATP-bindimJ domain-containinc protein 2 


EFN65530.I 


Camponolus floridanus 


2.00E-81 




R-SS III 14 


494 


ATP-dcDcndcnt C1d nrotcasc 


XP 002428757.1 


I'rJu ulm Immanm corooris 


2.00L-04 




R-SS Til 15 


699 


catheosin B endopentidase 


XP 002404475.1 


Ixodes scupuluns 


2. OUT. 32 




R 55, In 16 


186 


conserved hypothetical orotcin 


ABM53543.1 


uncultured beta aroteobacterlum CB\FD1 BAC clone 578 


i.oon-09 




R-SS 111 17 


818 


cyclic- AMP response element bindmp, protein isoform 2 


NP 001159940.1 


Nasonta vitnpenms 


2.00E-50 




R-SS TH 18 


■IIS 


DNA replication licensing factor MCM5 


P.FN70560 1 [Cammmntus fluiijunus. 


2.00E-44 




R-SS TH 19 


705 


DNA directed RNA polymerases 1, II, and III 1 7. 1 kDa polypeptide 


XP 002125189.1 


Pcdiculus humanus corporis 


1.0OE-39 




R-SS TH 20 


849 


DnaJ-like oiotein subfamily C member 2 


EFN63844.I 


Camponolus floridanus 


8.00E-51 




R-SS III 21 


844 


Exostosin-3 


EFN78322.I 


Harpesnalhos sallator 


2.001.-5 7 




R-SS TH 22 


584 


(1K2I721 


XP 002074462 1 


Dni\ii[>hilu nillisloni 


2 01)1 -OS 




R-SS TH 23 


■117 


(.jl.lOSSMb 


XP 002015686 1 


Drosophila pseudoobscura pscudoobscura 


2.00E-23 




R-SS TH 24 


398 


GroEL 


ACC54436 1 


Ciindidulux Pnrtieru tilrvnidiiliirum 


2 00F-50 




R SS TH 25 


245 


heat shock protein 90 


AUG03469.1 


Bemisia tabaci 


6.00k 34 




R-SS TH 26 


456 


hvpnlliclical protein AacL AAP1 O033S5 


XP 001656735 1 


Andes uexypti 


3.0OE-2I 




R-SS 111 27 


220 


hypothetical protein AcavDRAFT 4806 


LV 06212060.1 


Acidovvrax avenae subsp. avenae AI CC l')8bt) 


2.00E-21 




R-SS TH 78 


566 


hypnlhelical protein I Of 43 1729 


NP 001002182 1 


Damn rrri/i 


1 00F.-25 




R-SS TH 29 


779 


hypothetical protein TcasGA2 TC01 1939 


EFA09797.1 


Tribolium caslancum 


2.00E-27 




R-SS TH 30 


'S6 


hypothetical protein TcasGA2 TC030638 


EFA01434.1 


Tribolium casluneum 


1.00E-12 




R-SS TH 31 


482 


molecular chaperone DnoK 


YP 742734.1 


Alkalilimnicola ehrlichii MLHE-i 


8.00E-06 




R-SS TH 32 


Sll 


propionvl-CoA carboxylase alpha chain, mitochondrial-like 


XP 002937132.1 


Xenopus (Silurana) Iropicalis 


l.OOE-37 




R-SS TH 33 


311 


CG4364-PA 


XP 001600892.1 


Nasonia vitripennis 


I.OOE-33 




R-SS TH 34 


585 


ENSANGP00000013956 


XP 001604161 1 


Nasonia vitripennis 


4.00E-63 




R-SS TH 35 


435 


erect wine CG3 1 14-PF 


XP 001947009.1 


Acvrthosiphon pisum 


6.00F.-09 




R-SS TH 36 


312 


eukarvotic translation initiation factor 3 subunit 7 


XP 001950255.1 


Acvrthosiphon pisum 


4.00E-22 




R-SS TH 37 


368 


mitochondrial ribosomal orotcin L30 CG7038-PA 


XP 001942529.1 


Acvrlhosiphon pisum 


1.00E-09 




R-SS TH 38 


665 


Protein enabled 


XP 001943648 1 


Acvrthosiphon pisum 


2.00E-19 




R:SS TH 39 


. .470. .. 


signal recognition particle receptor alpha subunit (sr-alphal 


XP 001807129.1 


Tribolium castaneum 


3 00E-23 




R-SS I'll 40 


162 


SVH protein 


XP 001947154.1 


Acvrthosiphon pisum 


6.001. -09 




R-SS Til J! 


897 


tomosvnCG17762-PC 


XP 001952670.1 


Acvrthosiphon oisum 


2.00E-73 




R-SS III 42 


451 


vacuolar ATTase subunit C 


XP 001946227.1 


Acvrthosiphon pisum 


5.00E-59 




R-SS TH 43 

R-SS TH 44 


696 
876 


vasodilator-stimulated nhosnhonrotcin 

www ITttncrrinrflgi- 


XT 970236,1 

BAC57914.1 


Tribolium castaneum 

irirtnfcg;^ onmhina 


l,WE-18 

3.00E-17 




R-SS TH 45 


268 


SWl/SNF-related matnx-associated actin-dependent regulator of 


EFN84532.1 


Harpesnalhos sallator 


2.00E-32 




R-SS TH 47 


766 


Iransposase IS630 


CAL47051.1 


Listonella anxuillarum serovar 02 


8.00E-07 




R-SS TH 51 


1057 


unnamed protein product 


BAG64005.1 


Homo sapiens 


3.00E-71 




R-SS TH 52 


669 


vitellogenin-! 


BAA88077.I 


Plautla stall 


2 00F-09 




R-SS IH I 


570 


60S ribosomal protein L27a 


XP 001892460.1 


Brugia malavi 


7.00E-54 





* Genes randomly selected for the qRT-PCR validation study are in the shade; ** transcripts were putatively identified by 
homology from BLASTX search using a cutoff E-value less than IxlO 5 . 
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